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Abstract 

We formulate a dynamically triangulated model of three-dimensional Lorentzian 
quantum gravity whose spatial sections are flat two-tori. It is shown that the combina- 
torics involved in evaluating the one-step propagator (the transfer matrix) is that of a 
set of vicious walkers on a two-dimensional lattice with periodic boundary conditions 
and that the entropy of the model scales exponentially with the volume. We also give 
explicit expressions for the Teichmiiller parameters of the spatial slices in terms of the 
discrete parameters of the 3d triangulations, and reexpress the discretized action in 
terms of them. The relative simplicity and explicitness of this model make it ideally 
suited for an analytic study of the conformal-factor cancellation observed previously in 
Lorentzian dynamical triangulations and of its relation to alternative, reduced phase 
space quantizations of 3d gravity. 
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1 Motivation 



The approach of Lorentzian Dynamical TriangulationsQ (LDT) leads to a well-defined 
regularized path integral for 3d quantum gravity as was shown in || ||. The phase 
structure of this statistical model of causal random geometries has been investigated by 
Monte Carlo methods in the genus-zero case, where the two-dimensional spatial slices 
are spheres ||. Its maybe most striking feature is the emergence in the continuum 
limit of a well-defined ground state behaving macroscopically like a three-dimensional 
universe [7|, ||. This is in contrast with perturbative continuum arguments which 
suggest that in d > 3 Euclideanized gravitational path integrals are generically ill- 
defined because of a divergence due to the conformal mode. Since a Wick rotation 
from Lorentzian to Euclidean space-time geometries is part of the evaluation of the 
regularized state sums in LDT, one might expect to encounter a similar problem here, 
but this is not what happens. Instead, all indications point to a non-perturbative 
cancellation between the conformal term in the action (which still has the same struc- 
ture as in the continuum) and entropy contributions to the state sum (that is, "the 
measure"). It should also be emphasized that this cancellation is not achieved by any 
ad hoc manipulations of the path integral, for example, by isolating the conformal 
mode and Wick-rotating it in a non-standard way (in fact, it is quite impossible to 
isolate this mode in the non-perturbative setting of LDT). - Further discussions of the 
conformal-mode problem and its possible non-perturbative resolution can be found in 
ill- 

It is obviously of great interest to understand in a more explicit and analytic fashion 
how this cancellation occurs and how it gives rise to an effective Hamiltonian whose 
ground state is the one seen in the numerical simulations of 3d Lorentzian dynamical 
triangulations. 

Some progress in this direction has been made recently by mapping the three- 
dimensional LDT model to a two-dimensional Hermitian ABAB-matrix model [10, 11 1. 
This latter model has a second-order phase transition which is absent from the LDT 
model. This comes about because the matrix model naturally contains generalized 
geometric configurations which are not allowed in the original quantum gravity model, 
and which can be interpreted as wormhole geometries. The second-order transition is 
related to the abundance of such wormholes and the original LDT model corresponds 
to the weak-gravity phase of the matrix model below the critical value of Newton's 
constant .0 



1 See jl], ||] for recent reviews. 

2 The reason why a iwo-dimensional matrix model appears in the description of f/iree-dimensional quantum 
gravity is the fact that the 3d space-time geometries of the latter can be uniquely characterized by a sequence 
of 2d graphs representing the intersection patterns of the 3d triangulations at constant half-integer times 
t + 1/2. 
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Although the mapping to the matrix model has yielded some analytic informa- 
tion about the phase structure of 3d quantum gravity, the explicit transfer matrix has 
not yet been constructed. This is a desirable goal, because it would lead to a quan- 
tum Hamiltonian that among other things could be compared with already existing 
canonical quantizations of 3d gravity. Also, having a more detailed control over the 
combinatorics of the triangulated model would be extremely interesting in order to 
understand the precise cancellation mechanism between the conformal terms in the 
action and the entropic measure contributions. 

In the absence of a solution of the full three-dimensional model, one strategy is 
to formulate simplified models rich enough to capture the dynamics of 3d gravity but 
whose combinatorial properties at the same time are sufficiently simple to allow for 
an explicit solution. There are two types of restrictions one can naturally impose on 
discretized Lorentzian space-times. The first are restrictions on the allowed spatial 
geometries (these are two-dimensional simplicial manifolds built from equilateral Eu- 
clidean triangles) at integer proper times t, and the second are restrictions on the 
allowed three-geometries that interpolate between adjacent spatial slices at times t 
and t + 1. Note that the first type of restriction has a direct influence on the Hilbert 
space of the system, since the spatial geometries may be thought of as a basis in the 
position-representation (where "position" here stands for a spatial geometry). 

Such restrictions may or may not change the critical properties of the corresponding 
ensemble of random geometries. In a larger context, a task that still needs to be accom- 
plished is the classification of all possible three-dimensional LDT models according to 
their critical properties and dynamics as a function of the set of Lorentzian discretized 
space-times allowed in the state sum. (As already mentioned earlier, one may not only 
consider imposing restrictions on the class of geometries, but may also allow for gener- 
alizations.) In general, one would expect this large number of possible discrete models 
to fall into only a small number of universality classes. This turned out to be the case 
in two dimensions, where there are essentially two universality classes, depending on 
whether one allows space-time to grow so-called "baby-universes" or not |l2j (however, 



see [13] for some "exotic" variations). Of course, three-dimensional models of random 
geometries have been studied far less, and it is not a priori clear what structures one 
should expect to find. 

A number of "quantum-cosmological" LDT models of 2+1 gravity were considered 
in [14 1, see also the discussion in pi]. There, the number of three-geometries contribut- 
ing to the path integral was restricted by imposing symmetry constraints, reflecting 
homogeneity and isotropy properties of space. The space-time topology was fixed to 
I x T 2 , that is, with toroidal spatial slices. Imposing in addition (discrete) spatial 
translation invariance fixes the tori at integer-t to be locally flat and Euclidean. The 
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spatial geometries are then completely characterized by three numbers, namely, the 
torus volume and its (two real) Teichmiiller parameters. The reason for why one may 
still hope to capture the essential dynamical features of 3d quantum gravity this way 
- despite the drastic reduction in the degrees of freedom - is that canonical continuum 
considerations suggest that 3d quantum gravity has only a finite number of true phys- 
ical degrees of freedom (which are precisely the Teichmiiller parameters). Note also 
that the torus case is the simplest choice with non-trivial Teichmiiller parameters and 



also the one which has been most studied in the literature [15 



In the simplest and most restrictive model of such a torus universe one demands that 
also the spatial intersections at constant half-integer t should be lattice-translationally 
invariant jl4| (see Jlq] for related cosmological continuum models). This scenario is 
most easily implemented by choosing as fundamental building blocks tetrahedra and 
pyramids; see [17] for a generalization to 3+1 dimensions. Although it is straight- 
forward to work out the combinatorics of all possible interpolating three-geometries 
between two spatial slices, it turns out that the model does not possess an interesting 
continuum limit. This has to do with the fact that because of the strong symmetry 
restrictions there are too few "microstates" , that is, too few geometries contributing 
to the state sum at any given value of the action. More specifically, unlike in 3d LDT 
without such restrictions, the number of distinct triangulations of a given space-time 
slice At = 1 (corresponding to a single time step) grows only exponentially ~ e const - L 
with the linear size L of the spatial torus. Since this entropy term has to compete with 
the volume-suppressing exponentiated cosmological term from the Euclideanized ac- 
tion which is of the form e _AA * i2 = e _Ai2 , the state sum will always be dominated by 
geometries with effectively one-dimensional spatial slices as the number of tetrahedral 
building blocks goes to infinity. Any potential continuum limit is therefore unlikely 
to have anything to do with the original LDT model, and we must conclude that this 
cosmological model is simply not rich enough to study the conformal-factor cancella- 
tion and the effective quantum dynamics of three-dimensional quantum gravity. (For a 
general discussion of renormalization and continuum limits in dynamical triangulations 
approaches to gravity see reference |lj.) 

In the present piece of work, we will investigate an alternative, less restrictive cos- 



mological LDT model first introduced in [18]. Its spatial two-geometries are still given 
by flat tori, but the symmetry restrictions on the interpolating space-time geometries 
are relaxed. We describe this so-called hexagon model in Sec. ^. As usual, a given 
discretized space-time contributing to the propagator consists of a sequence of layers 
[t,t + 1]. For the hexagon model, the geometry of each such "sandwich" can be char- 
acterized as a tesselation of a regular 2d triangular lattice by coloured rhombi. The 
two-coloured graph dual to this tesselation is a superposition of two regular hexagonal 
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graphs describing the flat two-tori which form the space-like boundary of the sandwich. 
In Sec. H we compute the Lorentzian action for a sandwich geometry, together with its 
Euclidean counterpart. 

Our next task is the counting of all possible interpolating sandwich geometries for 
given torus boundaries. We show in Sec. ||| that the associated combinatorics is that 
of a set of "vicious walkers" on a 2d lattice with periodic boundary conditions. In the 
following section, we demonstrate that the contribution AS" to the action of a single 
sandwich is already essentially determined by the geometry of the flat tori which form 
its boundary. Still, there is a large number of microstates for given boundary data, 
and we prove that their number indeed grows to leading order exponentially with the 
torus volume, and not just linearly. In Sec. ^ we calculate explicitly the variables 
describing the flat spatial tori (for each torus, two real Teichmuller parameters and 
the two-volume) in terms of the data labelling a triangulated space-time sandwich. 
Since the latter are a set of discretized variables, it is of interest to see how they 
sample the usual continuous Teichmuller space of all flat tori. This is illustrated in 
Sec. by explicitly calculating the Teichmuller parameters for a set of geometries 
whose volume is smaller than a certain cutoff. We also include a sample plot of the 
associated moduli space, obtained by factoring out the large diffeomorphisms. Our 
conclusions are then presented in Sec. ^. - Appendix 1 contains details of the coordinate 
transformation between the discrete geometric parameters and the torus data, and 
Appendix 2 an asymptotic evaluation of the vicious-walker combinatorics relevant to 
the entropy estimate. 

2 The hexagon model 

The fundamental three-dimensional building blocks used in this model are (in the lan- 
guage of 0]) 3-1 tetrahedra, glued together pairwise, 1-3 tetrahedra, also glued together 
pairwise, and single 2-2 tetrahedra. The numbers i-j indicate that a tetrahedron shares 
i vertices with the spatial geometry at time t and j vertices with that at time t+1. As 
usual, the space-like edges of a tetrahedron all have squared length 1 (or a 2 in units 
of the lattice spacing a) and the time- like edges squared length —a (or -aa 2 ), where 
a > is real. The pairing of the 3-1 or 1-3 tetrahedra is obtained by gluing them along 
a (time-like) triangular face. 

The three types of building blocks are illustrated in Fig.[l]. Any three-dimensional 
"sandwich geometry" of height At = 1 we construct from these building blocks can 
be uniquely described by the intersection pattern that results when the tetrahedra are 
sliced in half at time t + 1/2 and the time-like triangles that are cut in the process are 
represented by one-dimensional links. We colour-code the links to distinguish whether 
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Figure 1: The three types of tetrahedral building blocks and their corresponding two- 
dimensional intersection patterns at t — 1/2. By definition of the model, both the 3-1 
and the 1-3 tetrahedra always occur in pairs. 



they come from triangles with tip at t (red links) or tip at t + 1 (blue links). Our 
three building blocks can thus be represented by blue and red double-triangles and by 
squares with alternating blue-red sides. Topologically, the intersection graph is again 
a torus. 



m=4 




l+w=12 

Figure 2: A regular triangular lattice serving as a "background geometry" at half-integer 
time t + 1/2. 

As already mentioned, we want to consider only amplitudes between flat two-tori. 
For the blue-red intersection picture this means that when the red links are simultane- 
ously shrunk to zero length, what must remain is a regular tiling of a torus by (blue) 
triangles, where exactly six triangles meet at each vertex (and similarly at time t + 1 
when we shrink away the blue links). A systematic way of generating intersection pat- 
terns with this property is as follows. Take a rectangular strip of a regular triangular 
lattice of discrete width I + w and height m, where the units are chosen in such a way 
that all vertices have integer coordinates^] (Fig.^). Let us for the moment take I + w 
3 The reason for splitting up the width into two integers will become clear in Sec. ||. 
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and m to be even, I = 21', w = 2w' , m = 2m', l',w',m' G Z + , since this will make 
it possible to identify the opposite sides of this strip without any twists to create a 
compact two-torus. (In a more general model, one may also allow for twists in either 
of these directions.) 

This regular lattice is to serve as a "background geometry" which we are going to tile 
with rhombic 2d building blocks so that no space is left blank. It is immediately clear 
that the intersections of the double-tetrahedra are blue and red rhombi. The squares 
that result from cutting the 2-2 tetrahedra will be "distorted" in this representation 
so that they can fit onto the triangular lattice. Since this can be done in two ways, 
we have a total of four rhombic tiles at our disposal (Fig.||a). It will be convenient for 
our purposes to adopt a dual notation where each rhombus is represented by a pair of 
crossing links. Each link connects two opposite edges of the rhombus and has the same 
colour, see Fig.|3|b. The rhombi can only be put onto the lattice if the colours of their 
edges (or their dual links) match pairwise at intersections. 



Figure 3: The four types of rhombic tiles (a), and their dual representation (b). 

The beautiful feature of this representation is the fact that any tiling of the strip 
(which of course must be compatible with its periodic identifications) automatically 
leads to in- and out-geometries which are flat, connected tori. The easiest way of 
seeing this is by following dual links (or pieces of dual links) of a given colour around a 
closed loop, where the loop must be such that no more lines of the same colour branch 
off into the loop's interior. Next, consider how this loop is represented in terms of 
triangles of the same colour which make up one of the adjacent spatial geometries. 
The pieces of straight coloured lines coming from the last two building blocks depicted 
in Fig]3|b do not correspond to any triangles at all. By contrast, the first two rhombic 
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Figure 4: These tiles appear as blue and red double-triangles in spatial slices of integer 
t. Their trivalent dual graphs form part of regular hexagonal graphs representing flat two- 
geometries. 

tiles correspond to a couple of adjacent triangles each in the relevant spatial geometry 
at integer t. Since the rhombic tiles can be put onto the triangular "background lattice" 
at half-integer t only with three possible orientations (see Fig.|5] below), it is clear that 
any of the loops introduced above corresponds to a sequence of exactly six triangles. 

In fact, these two types of rhombi can more properly be represented by dual trivalent 
graphs which are dual to the individual triangles, as shown in Fig.[|. It is then easy to 
see that each loop corresponds to a hexagon graph, with six "corners" of 120 degrees 
each, a property which has given the model its name. Translating this into a statement 
about the two-geometry, it means that there are six triangles meeting at each vertex, 
hence the geometry is everywhere flat. Thinking a little further along these lines, one 
can also convince oneself that it is not possible to obtain a flat triangulation of one 
colour that consists of two or more disconnected pieces. 



3 Gravitational action and transfer matrix 

The first step in constructing a path integral for our model is to determine which 
sandwich geometries At = 1 can occur and to compute their contribution to the action. 
Following ||, this action can be written as a function of the total numbers N31, N13 
and N22 of tetrahedral building blocks occurring in the slice, namely, 

S{At = l) = 

(N31 + N13) (irkyfa — 3k arcsinh — 3fc"v/a? arccos -r— — \/3a + 1 

v v3v4qT+1 4a + 1 12 

+N22(^k\/a + 2k arcsinh v^v 7 a + _ 4/ Cv /^ arccos V4a + 2 

V 4a + 1 4a + 1 12 

(1) 

where A and k denote the bare cosmological and inverse Newton's constants. The 
positive parameter a appearing in ([!]) describes the ratio between the squared lengths 
of the time-like and the space-like edges of the triangulation, Z^ ime = — a/g pace . Since 
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we will evaluate the state sums in the Euclidean sector of the theory, we need to Wick- 
rotate all of our Lorentzian discretized manifolds. As explained in detail in this is 
achieved by continuing a through the complex lower half-plane to negative real values. 
Let us for simplicity choose the standard value a = — 1, so that all edges have the same 
length. This gives rise to the Euclidean action 

S eu (Af = l) = (JV 31 +JV 13 )((-|7r + 6axccosi)fc+^A) 

+N 22 ((-2tt + 6 arccos -)k + ^| A ) 
= (JV31 +iV 13 )(-0.468fc + 0.118A) + #22(1.103*; + 0.118A) (2) 

Obviously in our model the numbers of building blocks of type (3,1) and (1,3) are 
always even. Note also that the action ^) contains boundary terms (the discrete 
analogues of the usual spatial integrals over the extrinsic curvature) in order to make 
it additive under gluing of subsequent layers of At = 1. The action without boundary 
contributions has the form 

Sgk(At = l) = (N 31 + N 13 )[(-tt + 3 arccos ~)fc+-±=A) 

+A r 22((-2vr + 4 arccos -)k + ^| A ) 
= (JV 3 i + JVi3)(0.551*! + 0.118A)+JV22(-1.359A: + 0.118A). (3) 

The partition function or propagator for a single time-step after the Wick rotation 
is given by 

G( 9l ,gr,At = l)= £ -U^^ (4) 

where the sum is over all possible sandwich geometries T interpolating between the 
two spatial boundary geometries g\ and g 2 , and C(T) is the order of the symmetry 
group of the triangulation T. As usual 0], expression (0) defines the transfer matrix 
T of the system with respect to the natural scalar product (51(52) = ctal)^ 1 ' 92 ' v ^ a ^ s 
matrix elements 

(g 2 \f\g 1 ):=G(g 1 ,g 2 ;At=l). (5) 



4 The combinatorics of the intersections 

In trying to characterize all possible intersection patterns that can occur (that is, 
all possible 3d sandwich geometries), it is convenient to break up the combinatorial 
problem into two steps. The first one is how to tile the triangular lattice with (identical) 
rhombi, and the second how to introduce a colouring on the resulting tilings (in the 
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form of drawing chains of dual coloured links onto the rhombi). 




Figure 5: There are three orientations in which a rhombus can be put onto the triangular 
background lattice. 

A rhombus can be put onto the lattice with three different orientations, which 
we call A, B and C. This is illustrated in Fig.[|, where the centres of the rhombi 
are indicated by small circles. Fig.|| shows a complete periodic tiling of a strip with 
(I + w,m) = (12, 4). (Opposite sides of this strip are to be identified.) What should be 
noted here is that the number of A-blocks as well as the combined number of B- and 
C-blocks per horizontal row is conserved as one advances in steps in the m-direction. 
Starting at some B- or C-block in row 1, one can therefore follow a "path" in vertical 
direction made up of some sequence of B- and C-rhombi until one reaches the upper 
end of the strip (shaded region in Fig||). 
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Figure 6: A tiling with A-, B- and C-rhombi (as usual, opposite sides of the strip are to be 
identified). The shaded region is a B-C-path of winding number 1 in the vertical direction. 

At this stage we will for simplicity impose a further restriction on the allowed 
patterns of rhombi, namely, that all B-C-paths should have winding number zero in 
the /-direction and winding number one in the m-direction. (This does not seem to 
impose serious restrictions on the in- and outgoing two-geometries, c.f Sees. [f| 0, but 
is a condition that could be relaxed, should this turn out to be convenient.) That is, 
following such a path between the lower and upper end of a strip, it should contain an 
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equal number of B- and C-rhombi, so that it closes on itself upon identification of the 
lower and upper ends of the strip. A configuration which violates this restriction but 
is nevertheless periodic is shown in Fig.[7| 




Figure 7: An example of a forbidden tiling. 

By shrinking all B-C-paths to zero width, one obtains a regular tiling of only A- 
rhombi (Fig.||, left). This reduced lattice can be thought of as a sublattice of the 
original strip in the sense that the chains of its dual A-links close onto themselves. If 
the number of B-C-paths was w/2, the reduced lattice will have width I and height 
m. Note also that these chains take the very simple form of straight diagonal left- or 
right-moving lines. A configuration of dual A-lines is most easily represented as a tilted 
square lattice (Fig.^[ right). The intersection points of the dual left- and right-moving 
lines will then have /-coordinate 0,2,4,6, ...in the odd rows and 1,3,5,... in the even rows. 




Figure 8: The regular tiling with A-rhombi obtained by deleting the B- and C-rhombi from 
Fig|| (left), and the corresponding dual tilted square lattice (right). 

As will become apparent in due course, working out the possible colour assignments 
for such dual A-lattices is an important part of the combinatorics of the model. They 
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are quite easy to enumerate. By assumption, any colouring of the dual chains has to 
respect periodicity in both I- and m-directions. The easiest way to obtain a consistent 
colouring is therefore as follows. Start at some (dual) vertex v at m = and colour 
the, say, left-moving A-line passing through v while following it around the lattice 
(keeping in mind the periodic identifications of the strip), until getting back to the 
original vertex. Using such a procedure, it is straightforward to see that the number of 
possible colourings for the entire configuration of A-lines depends on the integer d, the 
greatest common divisor of 1/2 and m/2. (For example, the lattice depicted in Fig. |8| 
has I = 6, m = 4 and therefore d = 1.) Picking an arbitrary vertex v at m = as origin, 
the colours of the left- and right-moving A-lines emanating from the first d vertices to 
its right (including v) may be chosen arbitrarily, with the remainder determined by 
periodicity. 

To count intersections of a certain type between the right- and left-moving A-lines 
(important for determining the action) it suffices to look at a fundamental diamond- 
shaped region oidxd vertices - this region will then be repeated lm/2d 2 times through- 
out the lattice. 




1 2 3 4 5 6 1 



Figure 9: The B-C-strands may be reintroduced by drawing paths onto the regular 
A-strip. In the figure, we have drawn the path corresponding to the shaded area in Fig 

We may now reintroduce the B-C-strands into this picture by drawing chains of 
links whose vertices in every row lie exactly in between the dual A- vertices. Starting 
from the initial row m = 0, we can again follow paths of dual B- or C-links by making 
at every vertex a choice of moving diagonally up to the left or to the right (Fig.^). 
More than one dual B-C-chain can pass through any one vertex, and neighbouring 
chains are allowed to share one or more links, but not to cross, so that their relative 
position along the ^-direction is preserved as we advance in m. (Obviously, to obtain 
configurations of the type depicted in Fig.|(| each of these paths must be "blown up" 
in the horizontal direction to width 2.) 

The combinatorics of the B-C-chains can be mapped onto a model of so-called 
vicious walkers (which are not allowed to touch) with fixed initial and final points, 
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by inserting columns of width 2 between every pair of adjacent B-C-chains^]. Various 
versions of vicious-walker models, differing in their boundary conditions for the paths 
and the underlying lattices, have been investigated in the literature. The most common 
choice is that of free boundary conditions in the /-direction, i.e. a lattice of effectively 
infinite width. The initial points of the walkers at m = are usually located at a 
minimal mutual distance A/ = 2 near the origin, for example, at /-coordinate 0, 2, 
4, ... , w — 2, and the final points at m max are either chosen freely or again grouped 
together at some point (lo, m max ) with /-coordinates Iq, Iq + 2, Iq + 4, ... , Iq + w — 2 



(see, for example, ]19| , |20| , 21] and references therein). An exception is the treatment 
by Forrester p2| . who uses periodic boundary conditions in the /-direction and walkers 
with equally (but not necessarily minimally) spaced initial and final positions. 

The case relevant for our 3d gravity model is that of periodic boundary conditions 
in both the /- and m-direction, and the combinatorial problem can be phrased as 
follows: Given three even integers I, m and w, how many ways are there of drawing 
w/2 indistinguishable vicious-walker paths with winding numbers (0, 1) onto a tilted 
square lattice of width I + w and height m ? The tilted square lattice is dual to similar 
lattices depicted in Figs.||, |9|, i.e. it has vertices on the horizontal axis at even parameter 
values 0, 2, 4, ... , l+w—2 and vertices on the vertical axis at values 0,2,4, ... ,m. Fig.|lC| 
shows a typical configuration of vicious walkers for (/, m, w) = (10, 10, 6). The vertical 
boundaries of this lattice are identified periodically as indicated in the figure, so that 
/ + w = 0. Each vicious- walker path starts on the lower horizontal axis at some point 
(1,0) and ends after m steps on the upper horizontal boundary at the point (l,m), i.e. 
at the point with the same horizontal coordinate. 




1 2 3 4 ... l+w 

Figure 10: A configuration of three vicious walkers, corresponding to (l,m,w) = (10, 10,6). 



Vicious walkers are imaginary creatures that will do vicious things to each other when meeting at a 
point and therefore avoid such encounters. 
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This situation can be viewed as a special case of that of a single random walker 
in an alcove of the affine Weyl group of type A w / 2 -i 23]. The reflections of the path 



at the walls of the Weyl chamber correspond in this case to the collision points of 
w/2 random walkers that move on a one-dimensional circle. The ensemble of walkers 
takes simultaneously steps of unit length along the circle, either to the right or the 
left. Mapping these onto diagonal upward steps on a tilted square lattice and requir- 
ing identical initial and final points for each walker on the circle leads exactly to a 
situation as depicted in Fig. [TO. Following pi 23], the number of non-intersecting 



path configurations for w/2 walkers with initial and endpoints A = (Ai, A2, . . . , X w /2), 
Aj 6 {0, 2, 4, . . . , / + w — 2}, ordered along the circle so that Ai < A2 < . . . < A to / 2 , is 
given by 



b(X,m,l,w)= J2 iH 
f,£t i= o 2X2 



m 

l+w 



(m/2) + ti + Xj — A 



(6) 



where t is a u;/2-tupel of integers. Note that because of the properties of the binomial 
coefficients, only a finite number of terms in the sum over t is non- vanishing. The 
presence of the determinants has to do with the fact that although all possible path 
configurations appear in @, the contributions from configurations with touching or 
intersecting walkers cancel appropriately by virtue of the alternating signs in the de- 
terminantal sum, in such a way that only the non-intersecting ones are left over. The 
largest term that can appear in any of the determinants on the right-hand side of (0) 
is always the product of the elements on the diagonal, namely, 



(7) 



corresponding to the independent product of w/2 periodic, free random walks consisting 
of m steps. 




5 Adding colour and estimating the entropy 

The hexagon model possesses a feature that will simplify the (asymptotic) analysis of 
the propagator. This has to do with the fact that each intersection pattern character- 
izing a sandwich geometry can be brought to a standard form by using a sequence of 
moves which affect neither the in- and outgoing 2d torus geometries nor the value of 
the action (Q). 

The basic idea is to move all B-C-paths to the far left of the strip. The elementary 
move necessary to achieve this is the flip of a left-right wedge (an adjacent pair of a 
dual B- and C-link) to a right-left wedge or vice versa (Fig.|ll]). Since any B-C-chain 
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Figure 11: The effect of a wedge flip on a B-C-path drawn on the dual A-lattice. 



has by assumption an equai number of B- and C-iinks, it can be moved to the left 
of the strip (where it assumes the form of a zigzag path) by a finite sequence of such 
moves. In the process, it will cross (pieces of) A-chains, but not other B-C-chains. A 
typical final result after applying this procedure to all B-C-chains of an intersection 
pattern is illustrated in FigjTJ. 



Figure 12: A rhombic tiling of standard form. 




T 
/ 



Figure 13: The rearrangement of the rhombi making up a fundamental hexagon region during 
a wedge flip and of the associated dual links or lines. 

In order to understand which consequences the wedge flip has for the geometry, let 
us analyze which action it translates to on the unreduced rhombic tiling. The region 
on the original triangular lattice which is affected by the wedge flip is confined to a 
set of six triangles forming a single hexagon, with six dual links emerging from the 
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sides of the hexagon. There are two ways of tiling this fundamental hexagon by three 
rhombi. In both cases, there is a piece of B-C chain entering at the bottom of the 
hexagon and coming out at the top, a piece of a right-moving A-chain entering at the 
bottom left and exiting at the top right side of the hexagon and a piece of left-moving 
A-chain entering at the bottom right and coming out at the top left. The wedge flip 
changes the tiling of the fundamental hexagon and the intersection patterns of the dual 
chains in the interior of the hexagon, without altering the dual links emanating from 
it (Fig©. 

To determine the effect on the geometry, one has to consider all possible colourings 
of the chains passing through the hexagon and how the wedge flip affects the blue and 
red dual link patterns. The first case is that where all lines have the same colour, say, 
blue, which implies that the hexagon represents a gluing of six 3-1 tetrahedra. The 
wedge flip only alters the way in which we mentally divide this set of six tetrahedra 
into three double-tetrahedra. Obviously this move does not in any way affect the 
three-geometry, and the geometries before and after the move should not be counted 
as distinct. However, for all other colour choices (there is a total of six) the three- 
geometry genuinely changes. As one can easily verify, the wedge flip in all of these 
cases corresponds to pulling a piece of a red chain (without vertices) across a blue-blue 
intersection or the other way round. This means that the individual red and blue dual 
graphs are squeezed and stretched in the process, but remain otherwise completely 
unaffected, and hence will correspond to the same two-geometries. 

Since moreover the number of building blocks of a certain type (3-1, 1-3 or 2-2) does 
not change under a flip move, we have proved our original assertion that the wedge 
flip leaves all two- and three-volumes (and hence the action) as well as the Teichmuller 
parameters of the two-geometries invariant^] We recognize here a simplified feature 
of the hexagon model, compared with the most general dynamically triangulated 3d 
gravity model [Q, ||], even if we restricted its integer-i slices to be flat tori. Namely, 
although the toroidal two-geometries forming the spatial boundaries of a space-time 
sandwich At = 1 by no means fix the three-geometry inbetween the two slices, they 
determine essentially uniquely the value of the sandwich action S(At = 1). In other 
words, there is a large number of interpolating space-time geometries for given, fixed 
boundaries, but they all contribute with the same weight e lS . A main task in solving the 
model is therefore the computation of the number of distinct interpolating 3- geometries 
between two adjacent flat two-tori. 

Before we can give a precise definition of this combinatorial problem, we still need 
to specify how we are going to parametrize the colouring of the rhombic intersection 

The wedge flip is an obvious candidate for a Monte-Carlo move in numerical simulations of the hexagon 
model; it will have to be augmented by moves that can change the r's and other physical variables. 
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pattern. Since we have already shown that any intersection pattern characterizing a 
three-geometry can be brought to standard form without affecting the individual red 
and blue torus geometries, we may without loss of generality think of the latter as 
subgraphs of this standard form. Using the notation introduced above, we will label 
the uncoloured standard form by three even integers (w, I, to), where w/2 is the number 
of zigzag B-C-chains on the left (whose individual height is to) and (I, m) is the size of 
the regular lattice of A-rhombi on the right. 

Let us now introduce a colouring by drawing closed dual lines onto this standard 
form, thereby producing a coloured standard form S. There are obviously w/2 inde- 
pendent vertical lines that can be drawn onto the B-C-columns. We will split them 
into w\/2 blue lines and W2/2 = (u> — w\)/2 red lines. It is clear that the order in 
which we colour these strands will affect the three-geometry, but not the individual 
two-geometries or the action. This results in a multiplicity for given in- and 

out-states, counting the number of possible orderings of blue and red vertical dual 
lines. 

We turn next to the colouring of the remaining dual lines, i.e. those that traverse 
the B-C-chains horizontally and the A-rhombi diagonally. The choice is restricted by 
the fact that the number of such lines which are closed (and therefore can be coloured 
independently) is exactly d for the right-moving and d for the left-moving lines. We 
will denote the numbers of blue right-moving and left-moving dual lines by a r and ai, 
and those of the corresponding red lines by b r = d — a r and b\ = d — a/. As before, 
the relative ordering of the blue and red lines in either direction leads in general to 
different three-geometries, but leaves the two- geometries and the action unchanged, 
thus contributing a factor ( f^) to the number of interpolating states. 

Putting all of these observations together, we can now rewrite the one-step propa- 
gator (H|) in a more explicit form. Using the essentially unique association (51,52) 
(l,m,wi,W2,ai,a r ) (c.f. Sec. || and Appendix 1), G now takes the form 

Of*,*, At = l) = Y. r^M(S)C +W2 )( ")(") e- S "' 5 > (8) 
_ C (6 ) \ wi ' V a r ' V ai > 



The number M{S) counts the distinct strip configurations that can be obtained by 
applying elementary wedge flips to the coloured standard form S, uniquely described by 
the six parameters (l,m,wi,W2,ai,a r ), and the relative order O of their dual coloured 
lines. (We are regarding configurations that differ by overall translations in the I- and 
TO-direction as equivalent. At any rate, this choice does not affect the remainder of our 
discussion.) 

In view of the discussion in Sec. |], we are interested in the continuum behaviour 
of (||), and in particular how the entropy contributions compete with the kinemati- 
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cal action S eu to result in an effective (continuum) action. The combinatorial part 
M(S) of (||) is very similar to that of the uncoloured problem described in connection 
with Fig.[lO], but the colour-dependence does now enter in a slightly subtle way. This 
again has to do with the slight overcounting present in our model (the subdivision of 
fundamental hexagon regions of one colour) which implies that not all vicious-walker 
configurations will correspond to distinct three-geometries. How often this occurs de- 
pends on both the boundary geometries gi and the relative order O of coloured lines 
of an individual sandwich geometry. It is clear that the overcounting will be most pro- 
nounced when the intersection pattern has very many dual links of one colour and very 
few of the other, because this will result in many local fundamental hexagon regions 



of one colour which are insensitive to wedge flips, cf. Fig. 13. 

Let us proceed on the assumption that - at least to leading order - the scaling 
behaviour of the entropy will not be affected by this overcounting. This is in part 
justified by the numerical investigations of ||, where we found that in the continuum 
limit, neighbouring spatial slices are strongly coupled, in the sense of having a similar 
total volume. Under this assumption, we can drop the sum over O in (||), and substitute 
the combinatorial factor by M = M, where 

M(m, I, w) := &(A, m, I, w) (9) 
A 

is the sum over all ordered w/2-tupels of initial conditions A for a set of random walkers. 
We would like to establish the behaviour of M in the limit as (m,l,w) simultaneously 
become large^]. This is not completely straightforward, since according to @ each b is 
a sum of terms that can be both positive and negative. 

For the purposes of this paper we will concentrate on establishing the leading di- 
vergent behaviour of M; we hope to return to the full analytic solution of the model 
in the near future. Is it the case that the leading behaviour is given by exp(c spatial 
volume), for some positive constant c? As explained earlier, this is needed for obtaining 
truly extended geometries in the continuum limit and a necessary prerequisite for a 
conformal-factor cancellation. The continuum limit that has been considered in most 
of the literature on vicious walkers is that in which both the width and the length of 
the lattice become large, but not the number of walkers. In this case one typically finds 
for the number M of walker configurations 

M = 2 numbcr of walkers X number of steps x ( subleading terms ) _ ( 1Q ) 

This is precisely the type of scaling behaviour we are looking for, but not the physical 



(i 



A natural canonical scaling ansatz is m,l, w — > oo while sending the cutoff a — > 0, in such a way that 
the dimensionful length variables M := am, L := al and W := aw remain fixed and finite. 
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situation we are interested in, since the paths here are "diluted" (even though their 
initial and final points may lie close together). 

A scenario of more immediate interest to us is again that considered by Forrester we 



already cited earlier [22|, who treats the case of N equally spaced walkers on a lattice 
of width /i with initial conditions Aj = iv, i = 1, ■■■,N, so that v = fi/N is the average 
distance between the walkers' paths. However, it should be noted that his class of 
path configurations is larger than ours: although he also considers periodic boundary 
conditions in the m-direction, an individual walker's path is not required to close on 
itself after m steps, but may wind around the lattice several times before doing so. The 
resulting closed paths will in general have (I, m)-winding numbers different from (0,1) 
(an example is shown in Fig. which by definition we have excluded from our current 
model. The evaluation of the combinatorics of this more general case turns out to be 
easier because it involves determinants of cyclic matrices, which can be simplified. For 



even m, v and odd iV the number of vicious-walker configurations is given by [22] 



iV " 1 2 1/72-1 2tt r> 

Af F (N,m,v) = II - E 2 m cos m (-(|U&)). (11) 

p=0 6=0 



However, it is easy to see that the overcounting involved in (|11|) does not alter the 
leading behaviour for large N and m compared to our more restricted case. Namely, we 
can divide the path configurations counted by ( |TT| ) into different sectors, depending on 
the lateral displacement A = X final ~ ^initial of the initial and final point of each walker. 
We have A. max = ±m, with the largest contribution to N F coming from configurations 
with A = (corresponding exactly to the path configurations with winding number 
1 we are counting in the hexagon model). The asymptotics of (|iT[) is much easier to 
handle, since the right-hand side is a product of sums over positive terms. Assuming 
an appropriate monotonic behaviour of Af F , it follows from the results in Appendix 2 
that the leading divergence for N, m — ► oo and fixed v > 2 is given by 

Af F {N, m, v) ~ c(v) mN , 2 cos(^) < c{v) < 2. (12) 

Coming back to the combinatorics of the function M, note that for a given lattice 
width and number of walkers there are (^"^jj^J t erms contributing to the sum 
(^), corresponding to all possible initial conditions for the walkers. Since this number 
grows at most exponentially in length (as opposed to volume), the leading exponential 
behaviour will coincide with that of the largest term in the sum, which is roughly 
speaking that of the configuration(s) where the walkers are maximally distributed (i.e. 
spaced out in the m-direction) over the available space and therefore get least in the 
way of each other. We can identify their average distance with the spacing v appearing 
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in Forrester's formulas. (Strictly speaking, v is of course an (even) integer, but we do 
not expect this to play a role in the continuum limit.) We may therefore conclude that 
for an average spacing v = (I + w) /w, and in the limit as m,l,w — > oo, the number of 
vicious-walker configurations behaves to leading order like 

M(m,l,w) ~ c{v) mw/2 x (subleading terms), (13) 



where c(v) is a constant of order 2, which appeared in (12). This shows that the hexagon 
model of three-dimensional quantum gravity has indeed enough entropy, compared with 
more restricted cosmological models that were studied previously. 



6 Teichmiiller parameters for triangulated tori 

Our next task will be to relate the parameters (l,m,wi,W2,a,i,a r ), introduced in the 
previous section to label a standard coloured three-geometry, to the variables describing 
the two spatial boundaries of that geometry. In this way we will achieve a clean 
separation of the geometric in- and out-data g\ and gi appearing as arguments of the 
amplitude G, and make contact with the standard parametrization of two-dimensional 
flat tori. 

It is a well-known fact that in the continuum any flat torus geometry g can be 
characterized by three numbers (f,7~i,T2), where v is the two-volume (in our case pro- 
portional to the number of triangles that make up the torus) and the Tj are the two real 
Teichmiiller parameters ^\ We would like to compute the geometric data (v^ ,r[ , t% ), 
i = 1,2, for the individual two-tori from the parameters (l,m,wi,W2,ai,a r ). 

Recall that the numbers of blue, red and blue-red rhombi at time t + 1/2 is given 
by N31/2, N13/2 and N22 respectively. The first two numbers give us directly the 
discrete two-volumes of the spatial slices at times t and t + 1, namely, = N31 and 
v ( 2 ) = iV"i3. (In order to obtain the volume in terms of the lattice spacing a, one needs 
to multiply with the volume of a Euclidean triangle, which is a 2 V3/4.) Taking into 
account the considerations of Sec. |I| on the division of the A-lattice into fundamental 
regions, one derives for these numbers the expressions 

Ttx Ivn 

N31 = w l^( a r + a l) + -J2 a r a h 
TTl 1 772 

N13 = W2 — (b r + k) + -pb r bi, 

TTl I TTl 

N22 = ^( w ±( b r + bi) + w 2 {a r + ai)) + -^p( a rbi + a/6 r ). (14) 



7 We will keep our options open about whether we eventually want to use these or rather the so-called 
moduli parameters, which label equivalence classes of Teichmiiller parameters with respect to the action of 
the mapping class group. 
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In order to determine the Teichmiiller parameters for the tori, we will identify 
for each colour (blue or red) a pair of oriented, closed geodesies, i.e. straight lines 
of minimal length on the corresponding torus at time t or t + 1. (Note that these 
do not have to coincide with any triangle edges.) Each pair consists of a circle with 
winding number (1,0) and one with winding number (0,1) when viewed as curves 
on the intersection pattern at t + 1/2 (the two winding numbers refer to the hori- 
zontal /, w-direction and the vertical m-direction). The winding numbers of closed 
curves in the spatial two-geometries can therefore be thought of as "inherited" from 
the three-dimensional geometry. These curves are clearly unique up to two-dimensional 
translations (if we think of the fiat tori at integer time t as being rolled out into the 
two-dimensional plane). For a given colour, their lengths and their relative intersec- 
tion angle are translation- invariant, and in one-to-one correspondence with the torus 
parameters (v, t\ , t%) . 





^2 



t+1/2 



Figure 14: An intersection pattern at time t+1/2 gives rise to a blue triangulation at time 
t and a red triangulation at time t + 1. We also show the two vectors Vi and V2 representing 
the closed geodesies of winding number (1,0) and (0,1) for both of these geometries. 



Let us look at a simple example in order to illustrate this procedure. FigjFj shows 
on the left a coloured intersection pattern in standard form representing a sandwich 
geometry, and consisting of 12 rhombi. It gives rise to two spatial boundary geometries, 
both of them flat tori with four triangles each. They have been cut open and are 
represented by two parallelograms in the plane, where the small black and white arrows 
indicate how their opposite sides are to be identified pairwise. Each of these two- 
geometries has oriented geodesies with winding number (1,0) and (0,1) respectively. 
They can be represented by a pair of vectors (Vi,V2) drawn onto the parallelogram. 
For ease of representation, the vectors are depicted on the right with their correct 
lengths and orientations, and with a common origin. 

It remains to compute the lengths and angles of these vectors from the data 
(l,m,wi,W2,a r ,ai,b r ,bi) characterizing the coloured intersection pattern in standard 
form. Representing the blue closed geodesic with winding number (1, 0) at time t as a 
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two-dimensional vector v} and the blue closed geodesic with winding number (0, 1) 
as a second vector one finds in terms of the discrete units inherited from the 

(I, m)-coordinate system at t + 1/2 



(i) 



V. 



(1) 



Tfi Tfl 

2d {ai ~ ar) '2d {ai + ar 



(15) 



Translating this into absolute units, one finds 



Vi 



v 2 « 



(^j( a l + a r) + Y' "IJ^ ~ flr ' ) ) a ' 



73 (°i -a r ),— — (a; + flr))a, 



id 



(16) 



where a again denotes the lattice spacing. The two vectors , V 2 characterizing 



the red triangulation at time t + 1 are obtained by substituting a; , a r , u>i in (|lq) by 
bi,b r ,W2- Note that the relations ( |l~6| ) imply the inequalities 



vf?>o, vg>o, v; ( ?>|v^|, V?? > |V 2 ( ?|. 



(17) 



As is shown in Appendix 1, the coordinate transformation between the six variables 
(l,m,wi,W2,ai/d,a r /d) and the components of the vectors Vj , subject to the con- 
straints 



V 



1(2) 



/(I) u'(2) 



'(1) 



1,2 ~~ K l,2 J K 2,l — K 2,l J 



(18) 



is one-to-one. 




T+l 



Figure 15: Standard representation of a flat torus with Teichmiiller parameter r = T\ + ir 2 
as a parallelogram in the complex r-plane. 



Up to a common rotation and a global rescaling, each pair of vectors (y[ , V^) can 
be identified with the vectors spanning the parallelogram in the standard representation 
of a flat torus with normalized area in the complex -plane (r® = r}^ -HV2 ), Fig.|l^. 
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It is straightforward to compute the SO(2)-rotation 



cos </>W sin (f>^ 
- sin cos </>W 



(19) 



which aligns the vector V x with the rj -axis. One finds 

/ClI a i ^ i \ w l\ if 11 a Vol , . . , 

cos0^ = -^ T y(— (^ + a r ) + — ), sm0^ = -^jj — (a z -a r ), (20) 

where denotes the length of the vector V[ , 

i W := llvf'll = ^(a 1+ar ) + ^) J + (^(a,-« r )) l B . 



(21) 



The corresponding expressions for cost/>( 2 ), sin(/>( 2 ), L^ 2 ) are obtained by substituting 
(ojjOdIOi) h-> (bi,b r ,w 2 ). 

We still have to rescale all lengths by a factor l/L^ 1 ) so that the rotated vector 



assumes length 1. Applying then both the rotation (|l9|), (20) and this rescaling 



to the vector V' 2 , we can read off directly the dimensionless Teichmiiller parameters 
. Collecting the geometric data describing uniquely the blue torus at time t, we 
have finally 



,(1) 



Im mwi 



( im 

{ a l a r-j2 + ( a « + a r) 2d 



(1) m (, 2 2\ l . / 2 

ri = i^(i)j2l^- a r) 5 + ^-«r) T ja, 



(i) m / , , . ^V3wi\ a 2 (i) , 00 s 

2 = i^LWpi a ^~ + (a ' + ar) ^~) a = X(LWP W ■ (22) 

./(fc) 



The map between the independent vector components V i • (or, equivalently, the sand- 
wich variables (/, m, wi,W2, a-i/d, a r /d)) and the (v^ k \r^) is in general two-to-one (see 
Appendix 1). 

One may also wish to reexpress the action in terms of these torus parameters. 
Because the transformation to these parameters is not bijective, this can only be done 
modulo a sign ambiguity. One first writes the counting variables Nm as functions of 
the two-vectors V'( fc ) and then in turn expresses the latter as functions of the torus 
parameters. Details of these calculations can again be found in Appendix 1. The 
explicit form for the sandwich action S eu (S) = S eu (v^\ t^ 1 \ v @\ t^) is given by (||), 
with 



A31 = v 



(i) 
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N 13 = v {2) 

" " 2 (r « +T f)) 2 + (r «_ r f)) 2 l v 2(Tl +Tl H w r (D +W r (2)J 

- ^((r 1 (1) ) 2 -(r 1 (2) ) 2 + (r 2 (1) ) 2 -(rf) 2 )). (23) 

The function 

^(Vffvg-Vgvg^ (24) 

can be expressed in terms of the (v^ k \r^) only up to a sign, namely, as the positive 
and negative square-root of 

« 2 - ^(ri" + r?? + <r»> - rf )>) - + „<^) 2 . (25) 



To To 



'2 '2 



The resulting action has a feature that may at first seem puzzling. Let us think 
of the exponentiated action for fixed boundary geometries as a matrix element (c.f. 
equation (||)), 

e -S»(,W,rW, fl P),rP)) = . (^(2)^(2)1^(1)^(1)). (26) 

Usually this kinematical "transfer matrix" can be written as i = e~ ah = l — ah + 0(a 2 ), 
where h denotes the discrete, kinematical Hamilton operator of the system^], and where 
we have reintroduced the discrete unit At = a for a single time step. 

For such an expansion to exist, in the limit of small At the configuration variables 
at a neighbouring spatial slice should always be expressible as v(t + a) ~ v(t) + av, 
r(t + a) ~ r(t) + af. However, looking at the explicit formulas for and r} 2 ^ 
(representing the variable t\ at two adjacent spatial slices) as functions of a; and a r , 
one sees that they always have opposite signs. It is therefore kinematically impossible 
to set them equal unless they both vanish (or, equivalently, a; = a r ). This can be traced 
back to the fact that a natural "pairing" between two neighbouring spatial geometries 
occurs in our construction by putting their two graphs together in an intersection 



pattern. As can be seen from the elementary example depicted in Fig. 14, the natural 
"dual" of a particular blue graph is given not by the graph itself, but by its reflection, 
where the relative orientation of the two coordinate axes has been reversed. In terms 
of the Teichmiiller parameters (since by definition we always keep r 2 > 0) this amounts 
to a map R : T\ \—* — 7*1, with a lowest-order expansion t = R + 0(a) for the matrix t. 
From this point of view, a more natural elementary time unit in our model consists of 
two time steps, a feature that has also been observed in other discrete models of 2+1 
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This is the "kinematical", and not yet the full, "effective" Hamiltonian, because we have not included 
any entropy contributions in the matrix element. 
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gravity f2§]. 



7 A discrete sampling of the Teichmiiller and 
moduli spaces 

Although we have established above a description of the geometry of spatial slices 
in terms of the standard Teichmiiller parameters, it is clear that in our model the 
parametrization will not be a continuous one, since the rW according to (^) are 
given as functions of certain discrete (integer) parameters characterizing triangulated 
geometries. From a continuum point of view our simplicial discretization therefore 
provides a discrete sampling of the Teichmiiller space T (which topologically is an 
R 2 ). It is interesting to look in more detail at how this sampling works as function of 
a given cutoff on the volume of either the space-time sandwiches or the spatial slices, 
say. 

In addition, one might be interested in the configuration space obtained from Teich- 
miiller space by identifying points that differ by the action of "large diffeomorphisms" 
(the so-called moduli space), i.e. those two-dimensional transformations that do not 
lie in the connected component of the diffeomorphism group of the torus. There are 
different suggestions of how these transformations should be incorporated into the 
quantum theory of 3d gravity, whether as exact invariances, as symmetries with a 



unitary action, or possibly not at all (see [26] for further discussions). We will not take 



any particular viewpoint here, apart from remarking that quotienting out by large 
diffeomorphisms usually leads to additional complications in the quantum theory . 

In order to determine the modulus of a given Teichmiiller parameter r = t\ + 1T2 , 
we have to map it by a sequence of modular transformations 

r-T + 1, r-^- (27) 

T T 1 

(the generators of the mapping class group of the torus) to a point in a fundamental 
region C T, usually taken to be the "keyhole region" defined by — 1/2 < t\ < 1/2 and 

| r |2= T 2 +r 2> 1. 

To get a qualitative idea of the nature of the discrete sampling of these two spaces, 
we have set up a small program that generates all values of (v, r) that occur for sandwich 
geometries with parameters < l,m,w < 10 (and with the lattice cutoff a set to 1). 
This is easy to implement, but does not quite amount to a systematic bound on the 
total space-time volume (which is proportional to V = (I + w)m). The range of V is 
between 4 and 200. Geometries with one or two vanishing spatial volumes «^ were 
not allowed. 
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Figure 16: Sampling Teichmiiller space, for maximal spatial volumes v max = 4, 14 (top row) 
and 24,34 (bottom row). 



We show in Fig. 16 a sequence of samplings of the upper half of the complex r-plane, 
with dots indicating the r-values that occur for a spatial slice of torus volume (the 
number of triangles, which is always even) v < v max , where v max = 4,14,24,34. The 
origin of the Teichmiiller space is clearly an accumulation point, with sampled points 
spreading out more and more with growing torus volume. The points are arranged 
symmetrically about the imaginary axis and seem to lie along various well-defined 
curves and straight lines through the origin. Note that any given point may occur 
at more than one volume; the information about this multiplicity is not included in 
our plots. Fig.|l7] shows half of the keyhole region with all moduli that occur up to 
torus volumes 32 (again, the left half of the region can be obtained by reflection). 



Comparing with the number of different points appearing in Fig. 16, it is clear that the 
modular parameters are highly degenerate. - It is clearly of interest to study these 
distributions more systematically and for larger volumes, and compare them with the 
natural measures on Teichmiiller and moduli space, but this would lead us beyond the 
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scope of the present paper. 
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Figure 17: A sample of moduli parameters inside the right half of the keyhole region. 



8 Conclusions 

We have introduced a dynamically triangulated model of three-dimensional Lorentzian 
quantum gravity whose spatial slices at integer times are flat two-tori. As we have 
shown, this symmetry restriction simplifies an exact analysis considerably: the spatial 
slices (and therefore the associated states of the Hilbert space) are labelled by just 
three parameters - two Teichmiiller parameters and a global "conformal factor" - and 
the evaluation of the model's entropy is related to that of a set of vicious walkers, 
which is a rather well-studied combinatorial problem. 

It should be pointed out that our model is not a "minisuperspace model" in the 
usual sense of the word ( "reduce classically and then quantize" ) , since the spatial slices 
at non-integer t are in general not translationally invariant. Therefore, more geometric 
degrees of freedom contribute to the superposition of space-time histories in the path 
integral than is obvious from the integer-t slices. (In this sense it is related in spirit 
to recently studied cosmological models in canonical loop quantum gravity, where also 
part of the reduction occurs only at the quantum level, see, for example, f28|.) As 
we already explained in Sec. |], the conditions on the space-time geometries cannot be 
too stringent in dynamically triangulated formulations of quantum gravity, in order to 
have a sufficiently large entropy and a potentially interesting continuum limit^J 

Crucially, we could show that the hexagon model does have enough entropy in the 

sense that the number of triangulations of a sandwich At = 1 with given boundaries 

9 One could argue that this is actually desirable from a physical point of view. The more restrictive a 
mini- /midi-superspace model is, the less likely it is that its dynamics is representative of that of full gravity 
(see j27| for further discussion and an explicit example) 
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to leading order scales exponentially with the volume of the slice. We also exhibited 
explicitly how the discrete triangulation data can be translated into the more familiar 
parametrization of the flat two-tori in terms of Teichmiiller (or moduli) variables. 

All of these properties give rise to the hope that the hexagon model will provide 
a link between the LDT formulation of three-dimensional quantum gravity and alter- 
native, reduced phase space quantizations. On the one hand there is a good chance it 
will lie in the same universality class as the full dynamically triangulated model (after 
all, the degrees of freedom we got rid of by working with flat spatial slices are those of 
the local conformal mode, which are known to be unphysical), and on the other hand 
the parametrization of the Hilbert space and the quantum Hamiltonian will be close 
to that of the canonical formulations. 

What remains to be analyzed is the precise nature of the continuum limit of the 
hexagon model. What is the subleading asymptotic behaviour of the state sum? Can 
we reconfirm our earlier conclusion [6|, [K| that the gravitational constant is not 
renormalized, but merely sets an overall scale? How is the divergence coming from 
the (global) conformal mode in the action compensated by a corresponding term in 
the entropy? What is the functional form of the effective continuum Hamiltonian in 
terms of the Teichmiiller parameters? What is its groundstate and is it identical to the 
groundstate seen in the numerical simulations reported in ||, ^, ||]? 

One may wonder how we can hope to make much analytical progress in the inves- 
tigation of a £/iree-dimensional statistical model. Firstly, the answer is of course that 
pure gravity is a very special type of theory in three dimensions, which is known to 
possess only a finite number of global metric degrees of freedom. Secondly, as we have 
shown, an essential part of the combinatorics of the hexagon model's transfer matrix is 
that of a iwo-dimensional problem of vicious walkers, about which a number of analytic 
results are already available. In summary, we think that a further investigation of this 
model is a promising avenue to pursue, both to advance our understanding of dynam- 
ically triangulated models and their continuum limits, and to achieve some degree of 
unification among the existing and rather disparate approaches to three-dimensional 
quantum gravity. 
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9 Appendix 1 



In this appendix we give some details of the coordinate transformations used in the 
main text in Sec.^. As a first step, we give the transformation law from (an independent 
subset of) the combinatorial parameters describing the discretized space-time slices to 
the dimensionful vectorial quantitites V . Taking into account the identites w\+w 2 = w, 
a i + h = d, a r + b r = d (with d := GCD(//2, m/2)), we choose as an independent set 
of the former the six variables (l,m,wi,w 2 ,ai/d,a r /d). Note also that not all of the 
components of the two- vectors V'W, i = 1,2, are independent, but we have the two 
constraints (|l8j). The set of transformation laws is given by 

* = l^O® + v£>), am = -L (vg> + yg>), 

9 V /(1) V' (1) 9 V ,(1) V /(2) 

,,/(l) 2 K 12 |/ 2 2 ,(2) 2 V 12 ^2,2 

awi = 2V 1 Y - , — , m , aw 2 = 2V 1 \ J - 



This is indeed a well-defined coordinate transformation, with Jacobian 

64 1 Vg } 



(28) 



2,1 



Next, let us establish the coordinate transformation that will enable us to rewrite 
part of the action in terms of two- volumes and Teichmuller parameters. Note first that 

a 2 ^iV 13 = V? 2) xV 2 (2) 

a^N 22 = vfxVf +V| (2) xVf ) . (30) 

These relations can be verified using the transformation laws ( p8|) from the independent 
variables (l,m,wi,w 2 ,ai/d,a r /d) to the vectorial quantities V. 

It is more involved to express the components of the two-vectors Vy* in terms of 
the torus parameters, subject to the constraints (|l8|). We need to invert the relations 

(o _ vg • vy } 

liv^ll 2 ' 
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T. 



(o _ x vf 

||V( W || 2 



„(0 = xV f * (31) 

This gives rise to certain algebraic expressions for the independent vector components 
(Vli i 2 j V'21 ' ^22 > ^li^ ' ) m t erms of the Teichmiiller parameters and the two- 
volumes, which we do not bother to write down explicitly here. The Jacobian of this 
transformation is 

J = 4ft(V; (1) x V? >)(vf } x Vf ^HV^lMlVpir 4 , (32) 



where f2 was already defined earlier in (|24|). Substituting the vector components into 
the equations (^), one obtains the expressions ( |23"|) given in the main text. The fact 
that this coordinate transformation is not bijective finds its expression in the fact that 
there are two regions in "V-space" where J has opposite signs, and which are separated 
by the hypermanifold defined by f2 = 0. In geometric terms, this comes about because 
for two pairs of two- vectors {V® ,V% }, i = 1,2 (where for each i the lengths and 
relative angle of v[ and V2 are uniquely specified by (t[ , t% , «W) according to 
(31)) there are in general two ways in which these two vectors pairs can be arranged 
relative to each other (for example, by specifying the angle between and ) which 
satisfy the constraints ( |l8| ) on the vector components, as well as the inequalities (|T7|), 
Note that this does not imply that if one of these V-configurations can be realized 
as a three-dimensional triangulation (i.e. corresponds to suitable discrete values of 
(l,m,wi,W2,ai/d,a r /d)) that the other one of the pair can too. In fact, this is in 
general not the case. 
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In this appendix we determine the leading asymptotic behaviour of the combinatorics 
of a set of equally spaced vicious walkers on a lattice with torus topology, as quoted 
in Sec. |5| of the main text. Following [22], the number of configurations of N vicious 
walkers whose initial and final coordinates are equally spaced at a distance v and who 
walk for m steps is given by 

N-i 9 m+i "/2-1 2 
Af F (N,m,v)= J] J2 cos m (-(l + b)), (33) 

where m and v are assumed even and N odd. The width of the lattice is therefore 
\x = vN and its height is m. We are interested in the asymptotic behaviour of (|33| ) as 
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m and ./V become both large. Assuming that v > 2, let us for some given p consider the 
corresponding sum over b on the right-hand side. It is clear that for large m, this will 
be dominated by the term with the largest absolute value of the cosine. For p < N/2, 
this is the term where b is minimal, i.e. 6 = 0, and for p > N/2 the term where b is 
maximal. We thus arrive at 



N-l 
2 



JV-1 



m f ~ 2^n- s m (2^) n cosm ( 2 ^+^)) 

p=0 



-,mN 



N-l 

2 



p=0 



uN- 



(34) 



where we have already dropped terms that do not scale exponentially with the volume. 
Next, let us make an estimate of the product in the square bracket for large N . Since 
for all N the argument of the cosine lies in the interval [0,ir/u] and since for large N 
we have 



,JV 



' 7T 

cos" I —(1 



N 



the leading asymptotic behaviour of N F can be estimated by 

(2 cos -) mN < M F (N, m, v) < 2 mN , v > 2. 



(35) 



(36) 



The argument leading to ( pq ) is not yet water-tight, since we first let m, and only 
afterwards TV become large. That there are no further terms contributing at the same 
order can be seen by establishing an upper bound for the potentially dangerous part 
of J\f F , 



JV-1 
2 



n 

p=l L 6=0 



Vt2 - 1 27T V ^ 

E — (f(| + 6)) 



N-l 

2 

<n 

P =i 



cos-(^)( 1 + (|-l 



cos m (^(l 



p )} 



N 



N-l 

2 

<n 

p=i 



cos™(^) ri+^-i)^w±i)) 



/ ,V , , 27T 7T 



JV-1 



n 

Lp=l 



COS 



,2ttp, 
vN' 



2 m 



(37) 



Letting now m and simultaneously become large, one sees that the term in front of 
the square bracket in the last line of (|37| ) scales at most exponentially with m, and not 
with the volume mN. 

For v = 2, a direct evaluation of (1331) leads to 



M F {N,m,u) 



(38) 
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which is independent of N. This happens because at v = 2 the walkers are "densely 
packed" and can only move in unison, effectively behaving like a single random walker. 
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